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Abstract: We review some recent results on how 'PT-symmetry, that is a simultaneous 
time-reversal and parity transformation, can be used to construct new integrable mod- 
els. Some complex valued multi-particle systems, such as deformations of the Calogero- 
Moser-Sutherland models, are shown to arise naturally from real valued field equations of 
non-linear integrable systems. Deformations of complex non-linear integrable field equa- 
tions, some of them even allowing for compacton solutions, are also investigated. The 
integrabilty of various systems is established by means of the Painleve test. 



1. Introduction 

There are many examples of non-Hermitian integrable systems in the literature pre-dating 
the paper by Bender and Boettcher which gave rise to the recent wider interest in 
non-Hermitian Hamiltonian systems. A well studied class of quantum field theories is for 
instance affine Toda field theories (ATFT) 

1 2 ^ 

£ = -a^(/)9'^0 + -2-^nfcexp(/3afc-0), (1.1) 
^ P k=o 

involving i scalar fields (p. The are integers often called Kac labels and the for 
k = 1, . . . ,i are simple roots with ao being the negative of the highest root. When the 
coupling constant /? is taken to be purely imaginary these models have interesting and richer 
features than their real counterparts. The classical solitons were found Q to have real 
masses despite the fact that the model is a non-Hermitian Hamiltonian system. Unlike as 
for real coupling, the scattering of their fundamental particles allows for backscattering such 
that the associated Yang-Baxter equations give rise to solutions in terms of representations 
of quasi-triangle Hopf algebras (quantum groups). For the simplest example, the ^i-model, 
corresponding to complex Liouville theory, a rigorous proof for the reality of the spectrum 
was found by Faddeev and Tirkkonen by relating it to the Hermitian XXZ-quantum 
spin chain using Bethe ansatz techniques. 
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In addition, integrable quantum spin chains of non-Hermitian type have been inves- 
tigated in the past for instance by von Gehlen [0]. The Ising quantum spin chain in an 

(2) 

imaginary field corresponds in the continuous Umit to the Yang-Lee model {A2 -minimal 
ATFT) 

1 ^ 

n = -J2''" + ^^t<^hi+i^''^h A,KGM, (1.2) 

i=l 

and may be used to describe phase transitions. Here N denotes the length of the spin chain 
and the erf, af are the usual Pauli matrices describing spin 1/2 particles and acting on the 
site i in the state space of the form (C^)®^. 

It is by now well understood how to explain the spectral properties of such models 
by means of "PT-symmetry, i.e. symmetry of the Hamiltonian and the wavefunction with 
respect to a simultaneous parity transformation and time reversal, pseudo-Hermiticity or 
quasi-Hermiticity |5|, |6|, |^. In addition is also established how to formulate a consistent 
quantum mechanical description via the definition of a new metric, although this is only 
worked out in detail for very few solvable models, e.g. Q. Interesting questions regarding 
the uniqueness of the physical observables still need further investigations and remain unan- 
swered. Once of the new feature is that unlike as for Hermitian systems the Hamiltonian 
alone is no longer sufficient enough to define the set of observables . 

For classical systems, which are the main subject of this article, one may also use 
PT-symmetry to establish the reality of the energy llC 



E= n[uix)]dx = - n[ui-x)]dx = / n^[u{x)]dx = eK (1.3) 



Note that unlike as for the quantum case the reality can be established from the Hamilto- 
nian alone, albeit together with some appropriate boundary conditions. 



1.1 PT-guided deformations 

Let us now turn to the question of how to use the above mentioned arguments to construct 
new consistent models with real spectra. In principle we could use any of them, but clearly 
to exploit PT-symmetry is most transparent, especially for classical models as we indicated 
in ([1.3|). Keeping in mind that the effect of a "PT-transformation is VT : x — > —x,p p 
and i —i, we may deform any "PT-symmetric function in the following way 

f(x)^f[-i{ixy], f[-i{ixyp'-\ /(^) + /[(i^)^] + /[(,^)-p-], (1.4) 

while keeping its invariance. The deformation parameter e S M is choosen in such a way 
that the undeformed case is recovered for e = 1. The same principle may be applied to 
derivatives of PT-symmetric functions 

dxfix) ^ U,s := d,,Jix) = -mur, dl, := ^r'^,,,, (1.5) 
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such that 



:= dlJ = -ie{if,flfL, (1.6) 

Ix 



fxxx;e '■ — ^x,ef — ^^(^/x) 



fxxx 1) [ "^^^ 



(1.7) 



Jx \ Jx 

and even to supersymmetric derivative of a 7'T-symmetric functions 

D = 9d, + de^ D,:=ed,,, + dg, (1.8) 

with 9 being the usual anti-commuting superspace variable. Remarkably it can be shown 
that the latter deformation can be carried out without breaking the supersymmetry of the 
models [11]. We shall apply these deformations in section ^ 



2. Prom real fields to complex particle systems 

The above mentioned principle appears at times somewhat ad hoc and often the only 
motivation provided is that such models are likely to have real spectra. However, in the 
context of integrable systems some complex particle systems arise very naturally when 
taking systems for real valued fields as starting points. 

2.1 No restrictions, ^-soliton solution of the Benjamin-Ono equation 

Let us consider a field equation for a real valued field u{x,t) of the form 

u{x,t) = -Y( AgM. (2.1) 

Chen, Lee and Pereira showed thirty years ago [12, ^] that this Ansatz constitutes an 



f-soliton solution for the Benjamin-Ono equation [14| 



Ut + UUj; + XHUj;^ = 0, (2.2) 

with Hu{x) denoting the Hilbert transform Hu(x) = ^ J^^d ^z^^^^ provided the z^. in 
(|2.1j ) obey the complex A^-Calogero equation of motion 

\2 

4 = y Y^izj - Zfc)-3, Zk G C. (2.3) 

This is certainly the easiest example to demonstrate of how complex valued particle systems 
arise naturally from real valued fields. 

2.2 Restriction to a submanifold 

Obviously we do not expect the above procedure to produce complex valued particle sys- 
tems when starting with any type of field equation. Dropping for instance in equation ( |2.2| ) 
the Hilbert transform and considering therefore Burgers equation instead will not lead to 
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the desired result. However, we may consistently impose some additional constraints and 
make use of the following theorem found more than thirty years ago by Airault, McKean 
and Moser |15|: 

Given a Hamiltonian H{xi, . . . , x„, ii, . . . , x„) with flow 

Xi = dH/dxi and Xi = —dH/dxi i = l,...,n (2-4) 

and conserved charges Ij in involution with H, i.e. vanishing Poisson brackets {Ij,H} = 0. 
Then the locus of grad I = is invariant with regard to time evolution. Thus it is permitted 
to restrict the flow to that locus provided it is not empty. 

In fact, often there are no real solutions to grad 7 = and one is once again naturally 
led to consider complex particle systems. We consider the Boussinesq equation, that is a 
set of coupled KdV type equations, as an example 

vtt = a(f^)xx + ft-Vxxxx + Vxx a,b £R. (2.5) 
Then the real valued field 

e 

v{x,t) = X^{x - Zk{t))~'^, AgM (2.6) 
fc=i 

satisfies the Boussinesq equation ( p.5| ) if and only if 6 = 1/12, A = —a/2 and obeys the 
constraining equations 

4 = 2V(z,-Zfc)-3 ^ 4 = _^|^, (2.7) 



zl = l- Y,{zj - Zk)-^ ^ grad(/3 - h) = 0. (2.8) 

Here I3 = Yl^j=i[^j /'^~^J2k^j ^j{^j~Zk)'^] and Ii = X]j=i two conserved charges in the 

j4£-Calogero model. In principle it could be that there is no solution to these equations, 
meaning that the imposition of the additional constraint ( |2.8D , besides the equation of 
motion (|2.7|), will produce an empty locus. However, this is not the case and some genuine 



non-trivial solutions may be found. For n = 2 a solution was already reported in |15] 

Zi = K+ V(t + k)2 + 1/4, Z2 = K- V(t + K)2 + l/4 (2.9) 

such that the Boussinesq solution acquires the form 

vix t) - 2X±^^l±l±3l±ll±. (2 10) 

Note that v{x,t) is still a real solution. Without any complication we may change k and 
R to be purely imaginary in which case, and only in this case, ( |2.10| ) becomes a solution 
for the "PT-symmetric equation (p.5|) in the sense that VT : x — > —x, t — > — t and v ^ v. 



Different types of solutions and also for other values of n will be reported elsewhere [16| 
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3. PT-deformed particle systems 

Having presented some examples of how to obtain complex many particle systems in a 
very natural way from real valued field equations, it appears less ad hoc to start directly 
by deforming some integrable many-body problems according to the principles described 



in section LI, having in mind that there might exist a corresponding real valued field 
equation. 

3.1 Complex extended Calogero-Moser-Sutherland (CMS) models 

The simplest way to deform a given model is just by adding a term to it along the lines 



indicated in equation (1.4). For a many body-system this was first proposed for the A£ 



Calogero model in [17| 



o2 , ,2 J 



with g,g G M, q,p G M^"*"^. The Hamiltonian TIbk differs from the usual Calogero model 
by the last term. There are some immediate questions to be raised with regard to (p.lj). 
Is it possible to have a representation independent formulation for TCbk'^ May one use 
other algebras or Coxeter groups besides Ai and Bg? Is it possible to use non-rational 
potentials? Can one have more coupling constants? Are the extensions still integrable? 



These questions were answered in [18|, where it was noticed that one may generalize the 
Hamiltonian TCbk to 

with A being any root system and the new vector /i = 1/2 ^ gaf{a-q)a, with /(x) = 1/x 

and V{x) = /^(x). It is not so obvious, in fact no case independent proof is known, that 
one can further re-write the Hamiltonian such that it becomes the standard Hermitian 
Calogero Hamiltonian with shifted momenta 



'Hf. = liP + il^f + I (3-3) 
and re-defined coupling constant 



2'^ 2 ^ 



.2 I 9s+aigi aeAs 

Here A; and A^ refer to the root system of the long and short roots, respectively. 

Thus we trivially have = r]~^hcai'n with rj = e~'^''^. Integrability follows then 
immediately by acting adjointly with rj on the Calogero Lax pair Lcai = [-^Cai, -^Cai], such 
that the new pair is obtained by L^{p) = Lq^i(p + ifi) and = Mcai- An interesting 
statement is obtained by computing backwards and allowing in ( |3.3D any kind of Calogero- 
Moser-Sutherland potential, i.e. V{x) = 1/x^, V{x) = l/sinh?x or V{x) = l/siv?x 

^^^ = \p^ + IY1 9ly{<^ ■q) + ifi-p- (3.5) 
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By construction the Hamiltonian (|3.5| ) corresponds to an integrable model, but it turns 
out H^] that the relation //^ = a'^g^ XloeAs ^(^^ ' l) + (^fsf YlaeAi ^{o: ■ q) is only valid for 
rational potentials. Thus without the ij? -term only the deformed version of the Calogero 
model remains integrable and not its generalizations. 

3.2 Complex deformed Calogero-Moser-Sutherland models 

Having seen that merely adding terms to complex Hamiltonians leads to rather simple 
models, we comment on some of the other possibilities indicated in (|1.4| ), which were 
explored in One of the symmetries of the CMS-models is its invariance with respect 
to the entire Coxeter group W resulting from the fact that we sum over all roots and 
the property that Weyl reflections preserve inner products. Interpreting now each Weyl 
reflection as a parity transformation across a particular hyperplane, we may try to seek 
models which remain invariant with regard to the action of a across these hyperplanes 
deformed version of the Weyl group associated with some newly defined complex 

roots a 

T-CpTCMS = Y^Te X] ■ + 2 ^ 9a^{^ ■ 9), (3.6) 

where m,ga G M. We outline the main features of the construction of the complex root 
system a € M" zM" with the desired features. First recall that to each simple root Oi 
there is an associated Weyl reflections ai{x) = x — 2ai{x ■ ai)/[aj). The aim is then to 
construct a new complex root system A in one-to-one correspondence to the standard one 
A, which may be recovered in the limit e = e — 1 ^ 

limaj(e) = Oj for a.i{e) S A(e),ai G A. (3.7) 

We define a T^T-Weyl refiection as da- := ct^^T, where the time reversal T has the effect 
of a complex conjugation. We may then use the action on a generic complex root a to 
determine their form 

CTa, {aj{e)) = Ga^T {^eaj{e)) + a^^T {ilmaj{e)) (3.8) 
= (Reaj(e)) - 10^^ (Iniaj(e)) (3.9) 
= — Reaj(e) 

As a solution to these equations we find 

ReQ;j(e) = R{e)ai and 

lim Rie) = 1 and 

with Aj denoting fundamental roots and /tj G M. Concrete examples for some specific 
algebras are for instance the deformed roots for A2 

diai{e) = —R{e)ai ^iI{e)X2 =: — ai(e), (3-14) 
-cTi(T2i7iai(e) = R{e)a2 =F iI{e)Xi =: a2(e), (3.15) 



ilniQj(e) 



(3.10) 
(3.11) 



Iniaj(e) = I{e)'^Kj\j, 
lim /(e) = 0, 



(3.12) 
(3.13) 
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or those for G2 

ai(e) = R{e)ai ± iI(e)A2, (3.16) 
a2(e) = R{e)a2 T i3/(e)Ai. (3.17) 

Having assembled the mathematical tools, we may substitute the deformed roots into 
the model ( p. 61 ) and study its properties. In jl^ it was found that the A2 and G2 deformed 
Calogero models can still be solved by separation of variables, analogously to the unde- 
formed case. However, some of the physical properties change, most notably the energy 
spectrum is different. The reason for this difference is that some restrictions cease to exist. 
For instance the wavefunctions are now regularized, such that no restriction arises from 
demanding finiteness. The original energy spectrum E = 2\uj\ (2n + A + 1) becomes in the 



deformed case |19| 



E^^ = 2\uj\ [2n + 6{Kf + Kf + £) + 1] forn,^€No, (3.18) 

with = (1 lb y/1 + 45f^/;)/4 and s, I referring to coupling constants multiplying terms 
involving short and long roots, respectively. 

In section 2 our starting point were real fields and we naturally ended up with complex 
particle systems, whereas in this section we started directly from the latter. It remains 
to comment on how the relation may be established in the reverse procedure. Having 
constructed the deformed roots we may compute the dual canonical coordinates q from 

a-q = q-a, a, g G M, d, g G M iM, (3.19) 



and subsequently simply replace in ( |3.1D the Hamiltonian H-prcMsiP, by 'H-pTCUsiP, Q, «) 
The freedom in the choice of the functions i?(e), /(e) may then be used to satisfy the con- 
straint grad/ = 0. The dual canonical coordinates for A2 are for instance easily computed 
to 

qi = R{e)qi + il{e)/3{q2 - gs), 

q2 = R{e)q2 + il{e)/?,{q^ - qi), (3.20) 
qs = R{e)q3 + il{e)/3{qi - ^2)- 

At this point it might still not be possible to satisfy grad/ = 0. As the construction outlined 
is by no means unique one still has the additional freedom to employ an alternative one. 



These issues are further elaborated on in |16|. 



4. Complex field equations 

Naturally we may also start directly by considering complex field equations. 
4.1 PT-deformed field equations 

Taking the symmetries of the Korteweg deVries (KdV) equation into account one may 
deform the derivatives according to ( |l.4D either in the second or the third term 

ut - 6uu;j;.s + Uxxx;fi = 0, e,;UGM. (4.1) 
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The first possibility, i.e. ji = 1, was investigated in |2C], leading to a not Galilean invari- 
ant, non-Hamiltonian system with at least two conserved charges in form of infinite sums 
allowing for steady wave solutions. Shortly afterwards the second option was investigated 
in that is e = 1, giving rise to Galilean invariant Hamiltonian system with at least 
three simple charges and steady state solutions. The question of whether these systems 



are integrable was thereafter addressed in [22| by carrying out the Painleve test. 
4.2 The Painleve test 



Let briefly summarize the main steps of this analysis proposed originally in |21|. The 
starting point is a series expansion for the field 



oo 

u{x,t) = Y,Mx,t)(t>{x,tf+'^, (4.2) 

fc=0 



with a being the leading order singularity in the field equation and \k{x,t) and (j){x,t) 
some newly introduced fields. The substitution of this so-called Painleve expansion into 
the partial differential equation (PDE) under investigation leads to recurrence relations of 
the general form 

X > ^'xx ) • • • ), (4.3) 

with / and g being some functions depending on the individual system under consideration. 
Solving these equations recursively might then lead at some level, say k, to g = 0. For 
that level we may then compute the right hand side of ( |4.3| ) and find that either / 7^ or 
/ = 0. In the former case the Painleve test fails and the equation under investigation is 
not integrable, whereas in the latter case A^ is found to be a free parameter, a so-called 
resonance. In case the number of resonances equals the order of the PDE the Painleve 



test is passed, since in that scenario the expansion (4.2) has enough free parameters to 
accommodate all possible initial conditions. A slightly stronger statement is made when 
the series is also shown to converge. In that case one speaks of the Painleve property of 
the PDE, which is conjectured to be equivalent to integrability. 

4.3 Painleve test for Burgers equation 



Instead of deformed KdV-equation (|4.l| ) let comment first on a simpler T^T-symmetrically 
deformed model, i.e. Burgers equations 

Ut + uUx-^£ = ii^Uxx;ii with K,e, fj, & M, (4-4) 



for which the Painleve test was carried out in 1 22 ] . In there it was found that leading order 



singularities only cancel for a = —1 and e = when taken to be integers. Keeping the e 
generic and starting with the lowest order, the recurrence relations lead to the following 
equations 

order — (2e + l): Xq + i2eK(j)x = 0, 

order — 2e: (pt^e,! + ^i4'x ~ i'^^^'xx = 0) (4-5) 

order — (2e — 1): d^ifptS^^i + Xicl)^ — iKe(pxx) = 0. 
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This means that 

Xq = —i2eK(j)^, Xi = {ieKcp^^ — (j)^5s,i)/(pxj A2 = arbitrary, (4-6) 
such that we have already one of the desired free parameters. A more generic argument 



can be used |22] to derive a necessary condition for a resonance to exist 

il^'-^e'Xrir + l)(r - 2)k^02^ = 0. (4.7) 



The requirement is that the parameter becomes free, which for (4.7) is obviously the 
case when r = —1,2. The values 2 was already found in ( [4.(^ ) and —1 corresponds to the 
so-called fundamental resonance, which seems to be always present. Thus according to 



the strategy outlined in section 4.2, we have the desired amount of free parameters and 



conclude that the deformed Burgers equation (W^ with e = /i passes the Painleve test. 



For the case e = 2 the convergence of the series was proven in [22| and concluded that 
the system even possess the Painleve property and is therefore integrable. In that case the 
expansion becomes 

One may even find a closed solution when making the further assumption u{x, t) = C,[z) = 
C(x — vt) and thus reducing the PDE to an ODE 

= . ■'V3(2.„«)V3 i'_-4';fa) + ai'(x) ^ 

c Ai{x) + Bt{x) 

with X = e^'^^^{vz — c){2vhc)^'^^^ and Ai, Bi denoting Airy functions. 

For the deformed KdV-equation ( |4.1| ) is was concluded that they are only integrable 
when e = fi, albeit the Painleve series was found to be defective meaning that is does not 
contain enough resonances to match the order of the PDE. 

4.4 Compactons versus Solitons 

A further interesting "PT-symmetric deformation was proposed for the generalized 
KdV equations which are known to possess compacton solutions 

ni,m,p = -777^ - -^uPiin^r. (4.10) 
/(/ — 1) m — 1 

The corresponding equations of motions are 

Ut + u^-^ux + gi'^uP-^^^-^ [p{p - l)ut 

+2pmuUxUxx + m{m — 2)ti^u^^ + mu'^UxUxxx] = 0. (4.11) 



The system ( 4.1(]| ) is yet another generalization of the generalized KdV system and the 



■pT-symmetric deformation of the KdV equation suggested in [10|, i.e. /i = 1 in (fl.lj), 



corresponding to Ti.i2p and Ws^^+i^O) respectively. These models were found to admit 
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compactor! solutions, which, depending on l,p or m, are either unstable, stable with in- 
dependent width A and amplitude (3 or stable with amplitudes /5 depending on the width 
A. An interesting question to ask is whether these systems also allow for soliton solutions, 
which would be the case when they pass the Painleve test. The test was carried out in |25 
and the results are summarized in the following table: 





compactons 


solitons 


I = p + m 


stable, independent A, f3 


no 


2 < / < p + 3m 


stable, dependent A, (3 


yes 


I <2 or I > p + 3m 


unstable 


no 



Interestingly, it was found that there is no distinction between the generalized KdV equa- 
tions and ■pT-symmetric extensions of generalized KdV equations with regard to the 
Painleve test. 

5. Conclusions 

We have shown that PT-symmetry can be used as a guiding principle to define new inter- 
esting models, some of which are even integrable. Some complex particle systems, possibly 
restricted to some submanifold, are shown to be equivalent to some real valued fields obey- 
ing non-linear field equations. We also investigated some complex field equations, which 
for certain choices of the parameters involved turn out to be integrable admitting soliton 
as well as compacton solutions. 
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